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Abstract 

We study by molecular dynamics simulations the wetting/dewetting transition and the dependence of the 
free energy on distance between plates that contain both hydrophobic and hydrophilic particles. We show 
that dewetting and strength of hydrophobic interaction is very sensitive to the distribution of hydrophobic 
and hydrophilic domains. In particular, we find that plates characterized by a large domain of hydrophobic 
sites induce a dewetting transition and an attractive solvent-induced interaction. On the other hand, a ho- 
mogeneous distribution of the hydrophobic and hydrophiUc particles on the plates prevents the dewetting 
transition and produces a repulsive solvent-induced interaction. We also present results for a kind of "Janus 
interface" in which one plate consists of hydrophobic particles and the other of hydrophiUc particles show- 
ing that the inter-plate gap remains wet until steric constraints at small separations eject the water molecules. 
Our results indicate that the Cassie equation, for the contact angle of a heterogeneous plate, can not be used 
to predict the critical distance of dewetting. These results indicate that hydrophobic interactions between 
nanoscale surfaces with strong large length-scale hydrophobicity can be highly cooperative and thus they 
argue against additivity of the hydrophobic interactions between different surface domains in these cases. 
These findings are pertinent to certain protein-protein interactions where additivity is commonly assumed. 
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I. INTRODUCTION 



The hydrophobic interaction is one of the major driving forces in various self-assembly pro- 
cesses such as protein folding, the formation of membranes and micelles, molecular recogni- 
tion and surfactant aggregation. [TD l2l El SI [H The nature of the hydrophobic effect is length- 
scale dependent. While hydrophobic ity at small length- scales is associated with small distor- 
tions of the hydrogen bond connectivity between the water molecules, large-scale hydropho- 
bicity is driven by a significant disruption of the hydrogen bonds network of interfacial wa- 
ters. [|6l 111 m HI [121 HH [121 [121 [131 The crossover from small to large scale behavior oc- 
curs when the radius of the hydrophobe is about 1 nm. One of the signatures of large-scale 
hydrophobicity is the existence of a dewetting (drying) transition; on bringing two large hy- 
drophobic plates from far apart into contact, a critical distance exists below which the wa- 
ter confined between the two plates is unstable in its liquid state and evaporates. The dewet- 
ting transition has been predicted theoretically JH [U] [141 [13 [IH and observed in computer 
simulationsJIIl [ni[Il[Il[20l[nil2l|2llM|251[2i|271[2l|2^ In a recent experiment us- 
ing interfacial-force micro scopyfBO], a nanoscale bubble was directly observed between super- 
hydrophobic surfaces. Nevertheless, the drying transition is sensitive to the magnitude of the 
solute-solvent interaction. This was demonstrated by showing that a small change of the attraction 
between the wall of a nonpolar carbon nanotube and water can induce a wetting transition|'22l. In 
this case, the transition point is characterized by a two-state model; in one state the channel is filled 
with water and in the other it is empty [[211 . As the solute-solvent attractive interaction increases, 
the dewetting transition becomes less pronounced and eventually disappears. [[TTl [T4[ 1271 [3T|| In 
some cases the kinetic pathway for dewetting in the cavitation transition has been studied using 
the powerful transition path sampling methodfT9l and the string method [|26l . The hydrophobicity 
of a surface is related to the value of the contact angle of a water droplet on that surface. Based 
on Young's equation, a surface with a contact angle larger than 90° is defined as hydrophobic. It 
is interesting to point out that, a hydroxylated silica surface with its partial atomic charges scaled 
by a factor of less than 0.4 was found to be macroscopically hydrophobic [[32ll . In addition, it has 
been argued that the mechanism for attractive mean forces between the plates is very different 
depending on the nature of the solute- solvent interaction [|33l. 

Both morphology and structure of the interacting hydrophobic surfaces are important to the 
existence and kinetics of dewetting. ||32l [34l [351 |36l In general, the surface of proteins is non- 



2 



homogeneous with respect to hydrophobicity, roughness and topology. For example, simulations 
indicate that protein surfaces can be characterized by hydrophobic regions that are heterogeneously 
"small" or "large" based on the length-scale dependence of hydrophobic hydration. [[371| . Thus, it 
is very challenging to predict the existence of the drying transition and determine its role in protein 
folding from the properties of the protein interface alone. In a recent study, we found that a simple 
hydrophobic scoring function, based on aligned hydrophobic surface areas, is not sufficient to pre- 
dict whether the assembly of protein units will exhibit dewetting or not. Il24l Further improvement 
of the scoring function should, therefore, include contributions from factors which were not taken 
into account, such as surface roughness and heterogeneity. 

A macroscopic theory relates the critical distance of the dewetting transition to the value of the 
contact angle of a water droplet on the surface, [[TT| 

" (P - n) + 2jijR^ 

where is the liquid/vapor surface tension for water, 9^ is the contact angle of water on the 
hydrophobic plate, P is the pressure of the liquid (water), Py is the vapor pressure of water, and 
Rm is the radius of the plates. For plates of nanoscale size, the term in the denominator (P — P„) 
is much smaller than 2'jiy/Rm and can be ignored, yielding, 

Dc^-Rm cos 9c (2) 

Thus, a knowledge of the contact angle of a heterogeneous surface is enough to predict the critical 
distance for dewetting. The contact angle 6* on a heterogeneous solid surface can be approximately 
predicted by the Cassie equation [[38ll39l . For a two-component surface, 

cos 9 = fi cos 9i + /2 cos 6^2 (3) 

where, fi is the fractional area of the surface with a contact angle of 9i. Simulations of droplets 
on heterogeneous patterned surfaces indicate that the Cassie equation holds for domains that are 
sufficiently small relative to the droplet. On the other hand, when the size of the heterogeneous 
domains is much larger than that of the droplet, there is a breakdown of the Cassie equation. [|40l 
I4TI |42| In the current work, we are interested in cases in which the size of the domains is small 
compared with the size of a typical macroscopic droplet. For example: under what conditions 
is the Cassie equation obeyed?; and how does the distribution of hydrophobic and hydrophilic 
domains (in terms of size and shape) determine the averaged contact angle 9 for a given /j? 
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Recent studies on hydrophobic hydration and hydrophobic interactions between hybrid po- 
lar/nonpolar nanoassemblies [l43l l44l l45l show that the pattern of nonpolar-site distribution is very 
important for the dewetting transitions to occur [|44l . Hydrophilic borders surrounding a nanoscale 
hydrophobic patch reduce considerably the patch's ability to repel water from the first hydration 
layer. In fact, even a single hydrophilic site at the center of the surfaces prevents complete drying 
of the confined region [|45l . In addition, it was found that an increase in the pressure of the bulk 
water blurs the difference between interfacial water density next to hydrophilic and hydrophobic 
surfaces [|45ll . However, none of these studies quantitatively addressed the correlation between 
the distribution of hydrophobic sites and dewetting as well as the strength of hydrophobic inter- 
actions. The related problem of how water behaves next to single heterogeneous surfaces has 
been addressed in a recent study using coarse grained modeling of the interace between water and 
heterogeneous surfaces [f36ll . 

It is also of great interest to investigate the extreme case where one plate is highly hydrophobic 
and the other hydrophilic. In an intriguing experiment Granick and coworkers investigated just 
such a system. They studied a Janus interface in which a water slab is trapped between a hy- 
drophobic wall on one side and a hydrophilic wall on the other [|46l and found that it prevents any 
macroscopic drying or cavitation of the liquid, which in any case would be strongly affected even 
by relatively weak van der Waals forces. 

Much work has been done on the non-additivity of hydrophobic interactions through calcu- 
lations of PMF between simple small hydrophobic solutes in dilute solution. [|47l |48l |49l [501 
[5D[5l|5SIMl55l|5S[57l[58l|5S[60lIthas been suggested that long range hydrophobic in- 
teractions caused by the many-body character of the PMF ||6T| is relevant for the energetics 
within [[62l [63l l64l or between [62J distinct "hot regions" of interacting protein surfaces. This 
is very important for the prediction and modulation of protein-protein interactions. Although the 
non-additive effect was found to be insensitive to the strength of the solute-solvent van der Waals 
interaction [|50ll , it was shown to increase with the size of the hydrophobic solutes Il50l [591 . The 
neglect of the length-scale dependence of hydrophobic solvation results in the failure to predict 
cooperativity for three-body hydrophobic association in current surface area based nonpolar mod- 
els ll65l . Previous studies of non-additivity were performed for simple hydrophobic solutes like 
methane, where the sizes of the clusters that form are small. Since proteins surfaces are het- 
erogeneous with mixed "small" and "large" hydrophobic regions, it is of interest to study the 
non-additive effect in the association of nanoscale hydrophobic assemblies in a heterogeneous 
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context. 

Our aim in this paper is to determine how the distribution of the constituent hydrophobic and 
hydrophilic particles determines the hydrophobic interaction between such amphiphilic surfaces. 
We perform systematic molecular dynamics simulations of different surfaces characterized by the 
same size and same number of hydrophobic and hydrophilic particles, which are distributed in 
different patterns. Our results show that the behavior of water confined between two identical am- 
phiphilic plates greatly depends on the distribution of the hydrophobic and hydrophilic particles 
on the plates. This effect is manifested by the existence or absence of a dewetting transition as well 
as by attractive or repulsive solvent induced interactions between the plates. We propose a param- 
eter, the average number of hydrophobic nearest neighbors of a hydrophobic particle, to describe 
the degree of clustering of the hydrophobic particles on the surface, and show that the behavior 
of confined water (dewetting or cavitation) is correlated with the value of this parameter. We also 
observe that for drying to occur the minimum area for clustered hydrophobic particles on the sur- 
face is (1.04x 1.04) nm^ and we call this the critical area for the drying transition. These results 
demonstrate a strong cooperativity in the hydrophobic interaction between surface hydrophobic 
domains of varied length scale. In addition, our results show that it is not possible to predict the 
critical distance for dewetting between heterogeneous surfaces by using the contact angle obtained 
by the Cassie equation (Eq.Q) in Eq.Q. We also carried out simulations of water confined to a 
Janus interface (where one plate is hydrophobic and the other is hydrophilic). In these simulations 
we find that the gap between these plates remains wet if the partial charges on the hydrophilic 
plates are sufficiently large. However, for small partial charges a drying transition is observed. 
In this case, the critical distance for dewetting can be obtained from the average of the critical 
distance for the pure hydrophobic system (where the two plates are hydrophobic) and that for the 
pure hydrophilic system. 

II. METHODS 

We studied the thermodynamics of the association process of two identical large amphiphilic 
plates composed of hydrophobic and hydrophilic particles. Each surface is represented by, a 
single-layered plate of, 49 atoms arranged in a 7x7 square lattice with a bond length of 0.32 nm. 
Each plate has 25 hydrophobic atoms and 24 hydrophilic atoms. The Lennard- Jones (LJ) parame- 
ters for the interaction between plate atoms on different plates are apu = 0.40 nm and epit = 0.50 
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kJ/mol, values lying in the same range as in our previous work. [|29l[66l We represent hydrophilic 
atoms in two ways. In the first approach, the hydrophilic particles are represented by the same LJ 
parameters as for hydrophobic particles, however, they have nonzero partial charges. Since all the 
surfaces we generated are electrically neutral, the number of positive and negative particles is the 
same. In one set of simulations we studied hydrophilic particles with partial charges of Q=±0.4 
e, and in another set with partial charges of Q=±0.8 e. In the second approach, we represented 
hydrophilic particles with a large well-depth of the LJ potential {^^f =1.30-2.0 kJ/mol), which is 

considerably larger than that for the hydrophobic particles (Cpf/^ =0.5-1.0 kJ/mol). Nevertheless, 
the LJ diameter was taken to be the same as that for the hydrophobic particles, a^it = 0.40 nm. We 
investigated five different patterns of hydrophobic/hydrophilic particle distributions on the plates 
(pattern I-V) (see Figure [TJ a)). 

The two plates were solvated in 1 147 water molecules. We chose the SPC/E model [l67l of 
water and used combination rules (arithmetic average for a and geometric average for e) to calcu- 
late the water-plate interactions. Analysis of the water density profile next to the plates suggests 
that the strength of the attraction between the water molecules and these plates is very similar 
to that of water and a hydrocarbon monolayer described at atomic level[68J. In addition, recent 
simulations [l69ll investigating the contact angle of water on the most hydrophobic surface studied 
here, e^f"^ =0.5 kJ/mol, find its value to be 119°. This value is similar to the values obtained for 
self-assembled octadecanethiol monolayer on silver llTOl (117°), and on gold llTTI (105°) surfaces. 

During simulations, the positions of the plate atoms are held fixed, interactions between atoms 
on the same plate are excluded, and the orientation of the two plates with respect to each other is 
parallel and in-registry. The alignment of the surface patterns with respect to each other is such that 
hydrophobic particles on one plate are superimposed on hydrophobic particles on the other plate. 
For hydrophilic particles, positive charged particles on one plate are superimposed on negative 
charged particles on the other plate. 

We used the Molecular Dynamics (MD) package GROMACS version 3.1.4[|72l to perform the 
computer simulations, with a time step of 0.002 ps. The bond distances and angle of the water 
molecules were constrained using the SETTLE algorithm [|73ll . The system was maintained at a 
constant temperature of 300 K and pressure of 1.0 bar using the Berendsen thermostatflTH. The 
electrostatic forces were evaluated by the Particle-Mesh Ewald method (with grid spacing of 0.12 
nm and quadratic interpolation) and the LJ forces by a cutoff of 1.0 nm. 

The potential of mean force (PMF) between the two plates was computed from the mean force 
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acting on each of the plates [l75l 17611 . Then the mean force acting between the plates along their axis 
of separation was integrated as a function of the distance between the plates, d, to yield the free 
energy profile. As the PMF represents only relative values, it was shifted such that the free energy 
of the state at the largest separation corresponds to zero. For each pattern (I-V), we performed 46 
simulations with different values of plate separation, d, ranging from 0.36-2.0 nm. At each value 
of d, the system was equilibrated for 2.0 ns and data was collected for 4.0 ns. At points where 
the force converged slowly (around the wetting/drying transition), the data collection stage was 
extended for an additional 5.0 ns. In the analysis of the hydration of the plates at each d for each 
pattern, the error in the quantities obtained from the simulations was estimated using the block 
averaging method llTTl . 

ni. RESULTS 

The dependence of solvation on the distribution of hydrophobic particles 

The behavior of water molecules between two amphiphilic plates depends on the distance 
between the plates, the distribution of hydrophobic particles on the surface, and the nature and 
strength of the interaction between the particles (hydrophilic and hydrophobic) and water. ll3T1l 
First we study plates in which the hydrophilic sites are charged spheres with either \Q\ = 0.4 e or 
\Q\ = 0.8 e to represent weak and strong hydrophilic particles, respectively. 

Figure. [2] (a) shows the density of water between two plates for the five different distribution 
patterns (shown in Figure [T]^a)) as a function of the inter-plate distance when the magnitude of the 
partial charge of the hydrophilic particles is equal to 0.4 e. Dewetting transition was observed in 
the inter-plate region for all the patterns except for pattern V, which is characterized by a uniform 
distribution of hydrophobic and hydrophilic particles. An absence of a drying transition for a 
uniform distribution of polar and non-polar sites was reported by Koishi et al[|44l|. In pattern V, 
the sharp decrease of water density near d = 0.72 nm is due to steric effects, i.e. at d = 0.72 nm a 
layer of water cannot fit between the two plates. Note that d is the distance between the center of 
mass of the particles on each plate along the z-axis. Thus, the available space between the plates 
is, approximately, d — Gpn. Interestingly, pattern IV with charged particles arranged in the center of 
the plates can dewet when \Q\ = OAe. Dewetting in this pattern might be due to surface dipoles of 
small magnitude formed by charged particles in which the distance between neighboring positive 
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and negative charged particles is small. This is in agreement with the finding that a silica surface, 
with partial charges scaled by 0.4, is macroscopically hydrophobic. [|32l The critical distance, Dc, 
for dewetting for each pattern is listed on Table |I} 

In order to correlate the geometrical pattern of each plate with its critical distance, we propose 
a parameter, the average number of hydrophobic nearest neighbors (A^„„) to characterize each 
pattern. We define Nnn as follows: 

1 ^ 

i=l 

Where Ni is the number of nearest hydrophobic neighbors of hydrophobic particle i, which has 
the maximum value of 4, and M=25 is the number of hydrophobic particles. We assume that the 
closer the hydrophobic cluster is to the center of the plate the larger will be the critical distance 
Dc for dewetting, but due to the small size of the plates in our system, Dc might not be sensitive 
to the distance between a hydrophobic particle and the center of the plate. Based on Equj4| we 
find a strong linear correlation between Dc and Nnn with the correlation coefficient r of 0.999 
(see Fig. |2](b)). Thus, as a measure of the degree of clustering for the assembly of hydrophobic 
particles on the plate, Nnn is capable of discriminating between the different patterns with respect 
to the solvation between two plates. 

When the magnitude of the charge of the hydrophilic particle is Q = 0.8 e, a similar trend was 
found for the effect of the distributions of hydrophobic particles on the solvation of the inter-plate 
region as in the case for Q = 0.4 e. However, a dewetting transition is suppressed in most patterns. 
Water depletion was found in patterns I-III and the extent of depletion is linearly proportional to 
Nnn for each pattern at the same plate-plate separation. Patterns IV-V stay hydrated until water is 
expelled because of steric effects. Note that water molecules are trapped more tightly in pattern V 
than in pattern IV, leading to a higher density of confined water than in the former, which is also 
consistent with their relative values of Nnn (see Fig. [3]). For all the patterns water density does not 
decrease to zero at d < 0.72 nm as in the case ofQ = 0.4 e, that is due to the strong electrostatic 
attractive interaction between water and the charged particles which reduces the distance for steric 
expulsion. 

Although we did not observe strong dewetting transitions for most patterns when Q = 0.8 e, 
we did observe stepwise cavitations for pattern II and pattern III. These patterns are characterized 
by multiple hydrophobic domains of different size distributed on the plate surfaces. Pattern II, 
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which has the second largest Nnn, contains two hydrophobic domains arranged in 4x4 (1.36x 1.36 
nm^) and 3x3 (1.04x1.04 nm^) square lattice domains. During the association of two plates, 
a big cavity is first formed between the larger hydrophobic domains (see Fig. [3] (b) when d = 
0.84 nm). This is followed by the formation of a second cavity between the smaller hydrophobic 
domains at smaller plate-plate distance (see Fig. [3](b) when d = 0.68 nm). For d <0.68 nm 
the hydrophobic domains are dry, and the water molecules between the hydrophilic domains are 
finally squeezed out into the bulk. It is interesting to point out that we did not find any stepwise 
cavitation for this pattern for weak hydrophilic particles (Q = 0.4 e). Instead, a large cavity 
forms without preference for the particles' hydrophobicity, covering hydrophobic and hydrophilic 
regions at the same time (see Fig.[3](b) when d = 0.84 nm). The stepwise cavitation we observe 
for pattern II is also observed for pattern III when Q = 0.8 e. However, in this case the cavities 
appear at smaller plate-plate distances, which almost can not contain more than one layer of water, 
compared with that for pattern II (see Fig [3] (c)). Thus, strongly hydrophilic particles restrain 
cavities locally to regions between hydrophobic domains and the cavities occur one after another 
based on the size of hydrophobic domains when bringing two plates from far apart into contact. 
These observed stepwise cavitations might be related to the size dependence of the amplitude 
of interfacial capillary-wave fluctuations as well as the probability of tube formation bridging 
vapor-film interfaces involved in the dewetting dynamics in the confined region. [fTSll The smallest 
hydrophobic domain area capable of inducing a cavity between the two plates, is found to be 
1.04x1.04 nm^ (see Fig|3](d)), which is consistent with the first cavitation in pattern III (Fig|3] 
(c) when d = 0.68 nm). Comparing the distances at which each cavity is formed during the 
association of the two plates for pattern II and pattern III, we find that stepwise cavitation for each 
pattern occurs at inter-plate distances which are linearly correlated with Nnn- Thus, the larger the 
value of Nnn for the entire pattern, the larger will be the gap distances of the stepwise cavitation. 

In order to investigate the sensitivity of these results to the nature of the particle's hydrophilic- 
ity, we repeated these simulations, for the case where the hydrophilic particles were uncharged 
(neutral LJ particles), but with a larger attractive interaction, epu, than for hydrophobic particles. 
In order to probe which values of can be regarded as hydrophilic and which hydrophobic, we 
performed additional simulations where all of the particles were taken to be the same and found 
that the threshold value of below which a drying transition is observed is 1.0 kJ/mol. Thus, 
we define particles with epu < 1.0 kJ/mol as hydrophobic and particles with epu > 1.0 kJ/mol as 
hydrophilic. Returning to the studies on heterogeneous surfaces, we performed simulations with 
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^^pit ~ 2-0 kJ/mol for hydrophilic particles and e^^" = 0.75 kJ/mol for hydrophobic particles. 
Fig. [5] shows the results for the two extreme cases of the distribution of hydrophobic/hydrophilic 
particles. We find that when hydrophobic particles are placed at the center of the plates as in 
pattern I, a strong dewetting transition is observed in the inter-plate region (red); whereas the 
inter-plate region remains hydrated if the distribution is uniform as in pattern V (green). These 
results are qualitatively the same as the case where the hydrophilic particles are represented as 
charged particles. 

In pattern I-V, the alignment of the two plates with respect to each other is such that the (X, Y) 
coordinates of the hydrophobic and hydrophilic particles on one plate are the same as on the 
other plate. It is also interesting to examine the behavior of water between two different surfaces. 
We study an extreme case of two plates which form a Janus interface [|46l (see Fig. [TJb)). One 
plate (plate 1 in Fig. [TJb)) is hydrophilic with positive and negative charged particles uniformly 
distributed on an 8 x 8 square lattice with the same bond length as for pattern I-V, while the other 
plate (plate2) is purely hydrophobic. The same simulations were performed for this case as for 
each of the other patterns. The density of water in this Janus -interface is shown in Fig. ^a) with 
respect to different partial charges Q. The Janus interface dewets when Q = 0.4 e and its critical 
distance Dc for dewetting is about 1.18 nm, which is a little bit smaller than the one (Dc =1.28 
nm) for Q = 0.0 e in which case both plates are purely hydrophobic. However, when the charges 
are increased to Q = 0.8 e, no dewetting is found between the two plates. Fig. |4] (b) and (c) 
show the density of water along the Z-axis of the simulation box (which is perpendicular to the 
plate surfaces) for inter-plate distances d = 0.72 nm and d = 1.24 nm, respectively. Only water 
molecules inside a rectangular box along the Z-axis of the simulation box with |X| <= 1.1 nm 
and |F| <= 1.1 nm (in the XY plane) are considered in these plots. (The origin of the coordinate 
system is the mid-point of the straight line connecting the centers of mass of the two plates. The 
hydrophilic plate is placed at the negative values of the z-axis.) These distributions indicate that in 
the presence of one purely hydrophobic surface, the strong hydrophilic surface {Q = 0.8 e) attracts 
water into the gap even when it would be geometrically impossible for this gap to accommodate 
one layer of water for Q = 0.4 or 0.0 e (see the sharp peak in red around z = nm in Fig. |4];b)). 
When Q = OA e, the density of water near the outside surface (toward solvent) of the weak 
hydrophilic plate is similar to that near hydrophobic surface with Q = 0.0 e (see the peaks in black 
and green with similar magnitude near z = —0.7 nm for d = 0.72 nm and those near z = —1.0 
nm for d = 1.24 nm). However, the density of water in the gap region increases for Q = 0.4 
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e compared to Q = 0.0 e for the plate-plate distance of 1.24 nm. These results are consistent 
with the observation of dewetting in pattern IV corresponding to Q = 0.4 e with hydrophilic 
particles distributed at the centers of the plates and hydrophobic particles on their borders. It also 
indicates that the weak hydrophilic plate with Q = 0.4 e is macroscopically hydrophobic. Thus, 
the behavior of water in the Janus interface depends on the polarity of the hydrophilic surface, i.e. 
the strength of the interactions between hydrophilic particles and water. 

Based on the simple macroscopic theory (see Equ. [T]), the average of A7 (A7 = — 7/^, cos Oc) 
should be used in the calculation of Dc for the two different plates in the Janus interface. To 
predict Dc for the Janus case with Q = OA e based on this simple theory, we repeated simulations 
for two identical, purely hydrophilic plates with partial charge Q = 0.4 e (same as platel in 
Fig.[TJb)). A dewetting transition was observed between these two plates and the critical distance 
Dc for dewetting was found to be in the neighborhood of 1 . 10 nm. We then determined the critical 
distance for the Janus case with Q = OA e from the average of for the pure hydrophilic plates 
(Dc = 1.10 nm) with Q = 0.4 e and for the pure hydrophobic plates (Dc = 1.28 nm). The average 
is 1.19 nm, essentially equal to the value (Dc = 1.18 nm) obtained from simulation of the Janus 
interface, indicating that the critical distance of dewetting in Janus interfaces can be predicted from 
the given Dc for pure hydrophilic plates (same to platel in Fig.[T](b)) and for pure hydrophobic 
plates (same as plate2 in Fig. [T](b)). 

A. Potentials of mean force 

Figure |6] (a) and (b) show the water induced PMF of the amphiphilic plates for the different 
patterns (TV), as a function of the plate-plate distance. For Q = OA e (see Fig. |6] (a)), the 
water induced force between two plates is attractive for pattern I-IV, while repulsive for pattern 
V, which has the smallest value of A^„„ among all five patterns and is the only one that does not 
display a dewetting transition. The difference in the water induced PMF for bringing the two plates 
from far apart to contact {d=OAO nm) between the different patterns is very large; it is about 230 
kJ/mol between pattern I and pattern V. While patterns III and IV display free energy barriers for 
dewetting, patterns I and II do not. At small distances, the water induced PMF is less negative 
for pattern I than for pattern II, probably because the water molecules like to stay in the gap due 
to their strong electrostatic interaction with charged particles on the edges of plates in pattern I 
compared with that in pattern II. When the partial charge of the hydrophilic particles is increased 
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to Q = 0.8 e, the water induced PMF is positive for almost all the patterns except for pattern I 
and pattern II in which the water induced force is attractive in a small range of inter-plate distance 
with a minimum at about d = 0.64 nm (see Fig. |6](b)). This corresponds to the cavitation or 
partial dewetting when two plates approach toward each other. The water induced repulsive force 
in the gap region for pattern V is very large in comparison with that for the case of Q = OA e. 
This means that it is very hard to remove water molecules from the inter plate region with strong 
hydrophilic particles. 

The solvent induced free energy of interaction (or the solvent induced part of the PMF) between 
two plates in the Janus interface as a function of their separation is shown in Figure [6] (c) for 
different partial charges on the hydrophilic plate. For both partial charges {Q = 0.8 or 0.4 e), the 
water induced interactions between the plates are attractive even though the Janus interface with 
Q = 0.8 e does not exhibit dewetting while the interface with Q = 0.4 e does. Nevertheless, 
for Q = 0.8 e the magnitude of the attractive interaction, as well as the shape of the curve which 
exhibits a solvent separated minimum and a barrier to remove this solvent layer, corresponds to 
an absence of dewetting. The difference in the solvent induced PMF for the Janus case and for 
the case with purely hydrophobic plates (Q = 0.0 e) (which also shows dewetting) increases 
significantly as the partial charge on the hydrophilic plate increases. This difference can be as 
large as 240 kJ/mol when Q is 0.8 e relative to that for Q = 0.0 e. 

IV. DISCUSSION AND CONCLUSIONS 

Our previous studies of proteins [|24ll indicated that although large matched and connected 
hydrophobic areas are correlated with a dewetting transition between two domains or oligomers, 
they are not sufficient to predict it. In this work, we aimed to determine a relationship between the 
magnitude of surface hydrophobicity and the spatial distribution of hydrophobic and hydrophilic 
domains on the surface. We determined the potential of mean force and through it the strength 
of the solvent induced interaction between two parallel identical amphiphilic plates, which should 
be regarded as an idealized model that might serve as a metaphor for protein inter-domain, or 
inter-oligomer interactions. The incorporation of hydrophilic particles in our model system was 
performed to mimic the effect of charged and polar side chains on the properties of the interface. 
We represented hydrophilic particles in two ways. The first is as particles with non-zero partial 
charges (±Q), and the second is by LJ particles with a large well-depth (e), significantly larger than 
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that for the hydrophobic particles. Since the number of hydrophobic and hydrophilic particles is 
constant (and the alignment of the different type of particles on the two opposing surfaces is in- 
registry, see Section |II]), the differences in the behavior of water in the gap must arise from the 
different spatial distributions of the hydrophobic and hydrophilic particles on the plates (at the 
given strength of water-particle interactions). 

We examined five different hydrophobic/hydrophilic particle distributions on the amphiphilic 
plates. Our results show that there are qualitative and quantitative differences in the behavior of 
the water for these different distributions. This is manifested by the existence or absence of a 
dewetting transition and by attractive or repulsive solvent induced interactions. Since the exis- 
tence of a dewetting transition is sensitive to the strength of the solute- solvent attractions [|3TI . 
we also studied how weak and strong hydrophilic particles affect the solvation of the inter-plate 
region for different plate patterns. In the case of weak hydrophilic particles (Q=OA e), most pat- 
terns exhibit a dewetting transition and the observed critical distance for dewetting, D^., varies for 
different patterns. The simulation results show a linear correlation between Dc and a proposed 
order parameter that describes the spatial arrangement of the particles on the plates. This order 
parameter, the average number of hydrophobic nearest neighbors (Nnn), is a measure for the de- 
gree of cooperativity for an assembly of hydrophobic particles on a surface. For example, for a 
pattern with a large cluster of hydrophobic particles at the center of the plates (pattern I), Nnn 
and Dc are both larger than for any other pattern. On the other hand, for a pattern where the hy- 
drophobic/hydrophilic particles are uniformly distributed (pattern V), A^^^ and Dc are both found 
to be smaller than for any other pattern. In this case, no dewetting transition occurs. Our proposed 
parameter, Nnn, seems to be able to discriminate between the different patterns and correlates 
highly with the critical distance for dewetting. For amphiphilic plates with strongly hydrophilic 
particles {Q= 0.8 e), the dewetting transition is suppressed in most patterns. However, the extent 
of water depletion in each pattern at the same plate -plate distance is linearly proportional to Nnn- 
Thus, Nnn is a good estimation of overall surface hydrophobicity. In addition, stepwise cavita- 
tions are found in systems with large hydrophobic domains in regions defined by the clustered 
hydrophobic domains. We found that the minimum area of the hydrophobic domain necessary 
to induce an adjacent cavity is 1.04x 1.04 nm^. Of course, the differences in the hydration of 
the inter-plate region corresponding to different patterns might lessen or disappear entirely if the 
hydrophilicity is made stronger. We did not investigate how the dewetting transition responds to 
misaligning the plates, but we expect it to be sensitive to their relative orientations. However we 
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did investigate the behavior of water between two plates forming a Janus interface [|46l consisting 
of one hydrophobic and one hydrophilic plate and found that the result depends on the polarity 
of the hydrophilic plate with the critical distance for dewetting being inversely correlated with the 
strength of hydrophilicity of the plate. The drying transition disappears entirely when the charges 
on the hydrophilic plate are sufficiently large (Q = 0.8 e). Our observation is consistent with an 
intriguing experiment performed by Granick and coworkers who investigated the hydrophobicity 
of a Janus interface. [|46ll . They found that the hydrophobic surface prevents macroscopic drying or 
cavitation of the liquid. This allowed them to focus on more intrinsic local properties of interfacial 
water near extended hydrophobic and hydrophilic surfaces, and to compare and contrast water be- 
havior in the different regions. Shear deformations produced by moving the hydrophobic surface 
resulted in very large noisy fluctuations consistent with the picture of damped capillary waves at 
the hydrophobic surface arising from partial dewetting. Film-spanning fluctuations that might lead 
to macroscopic dewetting between hydrophobic surfaces were suppressed by pinning of water at 
the hydrophilic wall. Our simulations give evidence that the critical distance of dewetting for 
the Janus interface can be predicted based on the simple macroscopic theory according to which 
it is the the average of Dc's for two pure hydrophobic plates (same as the hydrophobic plate in 
the Janus interface) and two pure hydrophilic plates (same as the hydrophilic plate in the Janus 
interface). 

In addition we also investigated the effect of the different patterns on the strength of the inter- 
plate interaction. This was done by calculating the PMF between the plates. In analogy to solvent 
induced interactions between hydrophobic particles, we found that for a pattern with a large hy- 
drophobic cluster, the inter-plate water-induced interaction is attractive (qualitatively, similar to 
the solvent induced interactions between homogeneous hydrophobic particles). However, for a 
pattern where the hydrophobic/hydrophilic particles are uniformly distributed, the solvent induced 
interaction is repulsive. This effect is substantial; the difference in the free energy change for the 
association process between these two patterns can be as large as ~230 kJ/mol. Physically, this 
qualitative difference in the induced potential of mean force can be attributed to the fact that for 
one pattern, it is easier, (negative induced PMF), on average, to strip off a water molecule from the 
plate interface than from around another water molecule in the bulk, while for another pattern it 
is harder (positive induced PMF). The solvent induced interaction between two plates which form 
a Janus interface is attractive (but less so for larger partial charges) even when the partial charges 
on the hydrophilic plate are large enough to suppress dewetting between plates. The free energy 
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barrier for dewetting for the strong hydrophilic system (Q = 0.8 e) reflects the free energy cost of 
stripping out the layer of water bound to the strongly hydrophilic surface as the plates approach 
each other. The PMFs of the Janus interface which display dewetting (Q = OA e) is qualitatively 
similar to other systems that exhibit drying. However, for the Janus interface with Q = 0.8 e, 
where the drying transition is absent, the curve of the PMF is different in shape and displays a 
minimum and a barrier that are associated with a solvent-separated layer. 

It is known that in order to describe stable native structure of known protein folds in a united 
(amino-acid) residue description, pairwise additive interactions are insufficient [TTSl 17911 . It has 
been shown [|80l [811 IHllllSII that only by including information about the many-body interactions, 
can one predict protein collapse or folding, similar to what is found in experiments. This cooper- 
ativity can arise from many groups in the protein. In this paper we demonstrated the importance 
of cooperativity for hydrophobic particles. Since non-additivity is likely to arise in systems with a 
strong solvent-induced effect, it is possible that the many-body effect found in proteins originates 
from the hydrophobic side-chain residues. 

In summary, we used molecular dynamics to study the thermodynamics of water confined 
between two amphiphilic plates, and found that different distributions of hydrophobic and hy- 
drophilic particles on the plates give rise to qualitatively different large scale water structures 
and water induced plate-plate forces. The results are sensitive to how strongly hydrophobic and 
hydrophilic the particles are. Our qualitative conclusions seem to be insensitive to whether the 
hydrophilic plates consist of particles which have coulomb interactions through partial charges 
or consist of particles that have strong LJ attractions with the solvent. Since for all patterns in- 
volved in the study of identical plates the number of hydrophobic and hydrophilic plate particles 
is the same, the results point to the breakdown of the Cassie equation and demonstrate that the 
hydrophobic interactions are strongly cooperative. We also investigated plate-plate interactions 
and the interplate large scale water structure in the Janus interface between a hydrophobic plate 
and a hydrophilic plate and found that when one of the plates is sufficiently hydrophilic it pins 
water molecules, eliminating the fluctuations that lead to drying. 
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TABLE I: The critical distance, D^, of dewetting transition for each pattern and its corresponding average 
number of hydrophobic nearest neighbor, Nnn (see Equ.[4]). 



Pattern 






I 


3.20 


1.06 


(1.04-1.08) 


II 


2.88 


1.02 


(1.00-1.04) 


III 


2.24 


0.96 


(0.92-1.00) 


IV 


1.96 


0.92 


(0.88-0.96) 


V 


0.00 


0.68 


(0.64-0.72) 
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Figure Captions 

Fig. [1] (a) The hybrid hydrophobic/hydrophilic plates with five different patterns of particle dis- 
tribution, pattern I-V. (b) Janus faced plates: plate 1 and plate2. The sphere in cyan represents 
hydrophobic particle, blue represent positive charged particle, red represent negative charged par- 
ticle. 

Fig. |2} (a) The density of water in the inter-plate region for pattern I-V as a function of inter-plate 
distance when Q = 0.4 e for hydrophilic particles, (b) The critical distance of dewetting transition 
as a function of the average number of hydrophobic nearest neighbors for each pattern. The red 
line is a linear fitting to the curve with correlation coefficient r of 0.999. 

Fig. |3} (a) The density of water in the inter-plate region for pattern I-V as a function of inter-plate 
distance when Q = 0.8 e for hydrophilic particles, (b) The view of a slab of water between two 
plates for pattern II. The water molecules are superimposed for 150 frames and are classified by 
different color based on their location. Water between opposite hydrophobic domains are repre- 
sented in green; water between charged domains are in red and blue; water in the bulk are in silver, 
(c) the same as (b) except for pattern III. (d) the same as (b) except for the system of two plates, 
each of which is composed by 9 hydrophobic particles arranged in 3 x 3 square lattice with a bond 
length of 0.32 nm. 

Fig. |4} (a) The density of the water molecules confined between two plates which form a Janus 
interface, for different hydrophilicity (magnitude of the partial charge Q) of the hydrophilic plate 
(platel, see Fig. [T](b)). Red corresponds to partial charge of Q = 0.8 e, black for Q = 0.4 e, 
and green for Q = 0.0 e in which case platel is purely hydrophobic, (b) The density profile of 
the water molecules along z-axis of the simulation box when the inter-plate distance d is 0.72 nm. 
The origin z=0 nm is the middle point of two centers of mass of two plates, and the hydrophilic 
plate is located at the negative values of the 2;-axis. Only water molecules inside a rectangular box 
along the Z-axis of the simulation box with |X| <= 1.1 nm and |F| <= 1.1 nm (in the XY plane) 
are considered in these plots, (c) The same as (b) except for d=\.2A nm. 

Fig. |5} The density of water between two plates where particles are the same (black, blue) and 
where particles are distributed in pattern I (red) and in pattern V (green). Here hydrophilic particles 
are represented as neutral LJ particles with e = 2.0 kJ/mol and the same is applied to hydrophobic 
except e = 2.0 kJ/mol. 

Fig. |6j (a) The water induced PMF of two hybrid plates for pattern I-V as a function of inter-plate 
distance when Q = 0.4 e for hydrophilic particles, (b) The same as (a) except that the charge of 
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hydrophilic particles is Q = 0.8 e (c) The water induced PMF of two plates which form Janus 
interfaces as a function of inter-plate distance with respect to different partial charge Q of the 
charged particles on platel (see Fig. [T] (b)). Black is for partial charge Q = 0.4 e, red is for 
Q = 0.8 e and green is for Q = 0.0 e in which case platel is purely hydrophobic and identical to 
plate2. 
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